Observables of a test-mass along an inclined orbit in a post-Newtonian approximated 
Kerr spacetime to leading-order-quadratic-in-spin 
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The orbital motion is derived for a non-spinning test-mass in the relativistic, gravitational field 
of a rotationally deformed body not restricted to the equatorial plane or spherical orbit. The 
gravitational field of the central body is represented by the Kerr metric, expanded to second post- 
Newtonian order including the linear and quadratic spin terms. The orbital period, the periastron 
advance, and the precession of the orbital plane are derived with the aid of novel canonical variables 
and action-based methods. 
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I. INTRODUCTION 

The motion of spinless particles in the gravitational 
field of a rotating body is an interesting and involved 
problem if the motion is not restricted to the equato- 
rial plane of the central body. Ever since the discovery of 
quasi-periodic oscillations (QPOs) and horizontal-branch 
oscillations in X-ray sources such as neutron stars, micro- 
quasars or black holes, theorists have tried to explain the 
presence of these oscillations by identifying their frequen- 
cies with those of nodal precession frequency of 'blobs' 
in the accreting disks or Keplcrian frequencies or perias- 
tron advance frequencies of the accreting rings. Since the 
mechanism that leads to these oscillations are associated 
with processes in the accretion disk, earlier calculations 
have been restricted to either close-to-equatorial plane 
or spherical or nearly circular orbits. Markovic in [5] 
studies the effect of equations-of-state of neutron stars 
on various orbital frequencies restricted to either equa- 
torial plane or close to equatorial plane. Merloni et al 
in [19] calculate the Lense-Thirring precession frequency 
of test-masses restricted to spherical orbits (constant r) 
and identifies them with the frequencies of QPOs from 
black hole sources. Sibgatullin in [4] calculates the nodal 
and periastron-rotation frequencies for nearly circular, 
non-equatorial orbits. Comparison of the observed and 
theoretical frequencies, once the relation between them 
confirmed, also act as a tool to test the general theory of 
relativity in strongly gravitating regions. Exact analytic 
solutions for arbitrarily oriented orbits are known in case 
of rotating black holes [1-3] but the Newtonian analog 
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or the post-Newtonian (pN) expansion of the results, say 
periastron shift, is not straightforward, see [4, 5]. Exact 
expression for the periastron advance of a test-mass in 
the gravitational field of a non-rotating black hole are 
well known, e.g. see [6-9]. Slow- motion or weak-field ap- 
proximations (pN approximations) of highly analytical 
or involved calculations allows one to have a meaning- 
ful interpretation of the results. In the references e.g. 
[6], [10], and [11] the motion of spinning binary systems 
were studied analytically and periastron advances were 
calculated to higher pN orders. The role played by the 
Carter constant [12] in various problems of motion re- 
lated with rotating black-holes is a quite interesting sub- 
ject too, see [13, 14], where Carter-like constants were 
introduced. This procedure will also be necessary for 
our calculations, for which we use Hill-inspired canonical 
variables [15-17] being better suited for the problem in 
hand than the often applied Boyer-Lindquist coordinates 
[18]. 

We, for the first time, without restricting to any special 
case, calculate the orbital period, the intrinsic periastron 
advance, and the precession of the orbital plane of a test- 
mass in an inclined, generic orbit in a Kerr spacetime 
field approximated to second pN order for the monopolc 
interaction, next-to-leading pN order linear in the spin 
parameter a, and the more important leading pN order 
quadratic in a 2 between the spinless test-mass and its 
spinning massive companion. For a highly inclined orbit 
with high eccentricity one needs to carefully handle the 
influence of the Carter constant which is absent in the 
case of equatorial orbits. We approach this problem by 
the most transparent way possible, which is by starting 
from an invariant scalar action in canonical variables al- 
lowing for a very clear definition and understanding of 
the observables deriving from it. For our calculations we 
use SI units. Our pcrturbativc method for making cal- 
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culations feasible is the post-Newtonian approximation 
technique, which expands field equations and equations 
of motion (EOMs) in inverse powers of c, the speed of 
light. Hence, if a term carries the power 4rj with n G No 
it is referred to as being of the ^th-pN order or Tj-pN in 
short where OpN is the usual Newtonian order. 



II. THE UNDERLYING ACTION PRINCIPLE 

The motion of the test-mass is governed by the binding 
energy, Ef,, which is conserved, and defines a Hamiltonian 
for generating the equations of motion. The total energy 
E of the system is given by 



E = mc 



E h 



(2.1) 



with m being the rest-mass of the test-body. We calculate 
the binding energy to 0(c -5 ), 
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(2.2) 



with the dimensionless spin parameter within 
in the range a G [—1,1], L z is the z-component of m's 
angular momentum, C is the Carter-like constant, r and 
p r being the radial coordinate and its canonical conjugate 
momentum respectively. The Carter-like constant, C, is 
defined as 
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(2.3) 
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We only consider leading order terms in a 2 , therefore, we 
can drop the _E 2 -tcrms being of higher pN order or next- 
to-leading-order. Notice that a 2 is positive because Ef, 
is negative for bound orbits. This Carter-like constant 
appears naturally as a separation constant when using 
the Hamilton- Jacobi approach to the pN-expanded Kerr 
metric in Boycr-Lindquist coordinates. We use a coor- 
dinate system that has canonical Hill-inspired variables 
(p r , L, L z ;r,u,£l), with u being the true anomaly, i.e., 
the angle between the position vector r = 7'n and the 



direction of the ascending node N = (e z x L) /\e z x L|: 
cos u = n ■ N and sin u = n • W with W = L x N/i. 
is the angle of the ascending node as measured from the 
x-axis of the non- rotating orthonormal basis (e x , e y ,e z ): 
cos £1 = N • e x and sin Q = N • e y , measuring the pre- 
cession of the orbital plane about L z which is tilted from 
the equatorial plane with the inclination angle i, the an- 
gle between the background z-axis and the orbital angu- 
lar momentum vector L, so L z = Lcosi. The standard 
canonical Poisson bracket relations arc 



{r lPr } = {u,L} = {n i L z } = l 



(2.5) 



with all other brackets being zero. In contrast to the 
classic Hill variables (f, G, H; r, u, ft) with G = r x v (v 
being the coordinate velocity vector), H = G z and (u, Cl) 
being conjugate to G and H respectively, we work with 
the momentum and its derived quantities p r = p • r/r, 
and L = |L| with L = r x p, i.e. our (u, f2) differs from 
(u,Q). They are equal in the Newtonian case (per unit 
mass), see [20] for variable relations between these dif- 
ferent frames. Notice that our orbital plane is defined as 
being orthogonal to L and not to G, where G is often 
denoted by L^r apart from a mass parameter. The cor- 
responding action integral about a complete revolution 
from periastron to periastron reads 



S = S(E b ,L z ,C,P,®,U) 



= -E b P + L z § + <j) p r dr + <j> Ldu. 



(2.6) 



It depends on the three constants of motion, E b ,L z ,C, 
and three orbital-completed variables 



P 



d£! and U 
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P being the orbital period, i> the precession of the or- 
bital plane per revolution, and U the intrinsic periastron 
advance. The relation between the Hill-inspired variable 
u and the Boyer-Lindquist coordinate 9 [21] is 



sin u sm i = cos t 



(2.8) 



Solving Eq(2.2) for p r and Eq(2.3) for L allows 
the straightforward calculation of the action, Eq(2.6). 
Within the context of our leading order calculation, L 
in cosi can be replaced by C. The next step is to calcu- 
late the remaining integrals in the action (2.6). We start 
with 

<j> Ldu = <j> \/C 2 — a 2 a 2 sin 2 i sin 2 u du 

= <j> Cdu - <j> -^a 2 a 2 sin 2 i sin 2 udu + 0{a A ) (2.9) 
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Since the integrand is of leading-order in a 2 we integrate 
over a closed Newtonian orbit, i.e., we integrate from to 
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2ir and drop all next-to-lcading-ordcr (NLO) corrections. 
The action, adapted to our problem, reads 



71" -2„2 / - L ' 



S = -E b P + L 2 $ + <p Pr dr +CU- —c l 2 a 2 [ 1 - -£ 



C 2 
(2.10) 



(2.11) 



The action principle tells us 

dS _ OS _ dS _ 
dC ~ dL z ~ dE b ~ 

giving rise to formulas for the orbital elements, 



$ = I p r dr - a\a 2 (E b )^ , 



(2.13) 



d f , a 2 7rG 2 Af 2 m / L 2 \ 
P= dEj Prdr+ ^^ V-t> ] - (2 ' 14) 



The integral over p r can be evaluated perturbatively in 
inverse powers of c, where p r is evaluated by inverting 
Eq(2.2) to 0(c~ 5 ). To make the integration easy we 
choose the method of contour integration employed by 
Sommerfeld in [22], 

<j) p r dr = —2iri [Res r= o(p r ) + Res r=00 (p r )] (2.15) 



The result reads, also see [6, 23]: 
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The two periastron shifts, AU and A$, the intrinsic and 
the one related with the precession of the orbital plane, 



respectively, are then giving by 

AU = ±-(U-2n) 

3G 2 M 2 m 2 6aG 3 M 3 
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We also recover the usual pN result when the motion 
takes place in the equatorial plane. In this case, C = 
L = L z and our precession angle $ agrees to leading pN 
order with Eq. (43) in [5], where it is called the nodal 
precession rate. Furthermore, the leading-order linear- 
in-a terms of AU and A$ fully agree with the results in 
[6] for the case of an inclined orbital plane when referred 
to the test-mass limit. For the period we have 
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In the case of equatorial motion the resulting shift, A<fr, 
(for this degenerate case) is given by 
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This quantity has been calculated by other methods [11] 
whose result coincides with ours, also see [5] to leading- 
pN-order. Notice that one can replace the Carter-like 
constant C by L in the Eqs. (2.17)-(2.20) because the a 2 - 
correction in C is lifting the pN order of the corresponding 
terms beyond our considerations. 
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III. APPLICATION TO MEASURED DATA 

As an application, we calculate the orbital el- 
ements of star S2 (treated as test-mass) orbiting 
the radio source Sagittarius A*, an assumed super- 
massive black hole with the mass M = 4 • 10 6 M© 
[24] and a spin parameter a = 0.44 [25] . We take 
the data of S2 star from SIMBAD. The data are 



Data label 


Measured value 


Gravitat. constant 
Solar mass 
Scmimajor axis 
Astronomical Unit 
Eccentricity 
Inclination angle 
Speed of light 


G = 6.67384(80) • 10" 11 m 3 /(kg • s 2 ) 
M Q = 1.98855(25) • 10 30 kg 
980 AU 

149.60 ■ 10 9 km 
0.881 

i = -48.1° 

c = 299792458 m/s 



The analysis of the orbital elements is as follows. 
We first give the Newtonian value labelled with N, 
then the pN monopole corrections to it with the la- 
bel 0, followed by the spin-orbit corrections (lead- 
ing and, formally, next-to-leading order) labelled by 1 
and finally the spin-squared corrections labelled by 2. 



the error bar for the period P = (15.56 ± 0.35)yr. 
Though we have used Hill-inspired variables to calcu- 
late the orbital elements only to some pN order, it 
will be interesting to see if these canonical variables 
would be advantageous to calculate the orbital elements 
non-perturbatively starting from a non-pN-approximated 
Kerr spacetime field as was done in Boyer-Lindquist coor- 
dinates by [26, 27]. At least in the case of equatorial mo- 
tion, the orbital elements calculated in Hill-inspired vari- 
ables should match the orbital frequencies given therein. 
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P (S2) 


AU (S2) 


Afl (S2) 


P N = 15.3484 yr 
P = 0.0012 yr 
Pi = -2.5 • 10~ 9 yr 
P 2 = yr 


AC/ N = 
AU = 701.115" 
AUi = -5.537" 
AU 2 = 0.008" 


An N = o 

Atio = 
Afli = 2.761" 
Afl 2 = -0.008" 



As one sees, the spin effects are negligible and within 



[1] E. Hackmann, Geodesic equations in black hole 
space-times with cosmological constant. PhD thesis, 
Bremen, 2010. 

[2] E. Hackmann et ai, "Analytical solution of the geodesic 

equation in kerr-(anti-) de sitter space-times," Phys. 

Rev. D 81 (2010) 044020. 
[3] E. Hackmann and C. Lammerzahl, "Observables for 

bound orbital motion in axially symmetric 

space-times," Phys. Rev. D 85 (2012) 044049. 
[4] N. Sibgatullin, "Nodal and Periastron Precession of 

Inclined Orbits in the Field of a Rotating Black Hole," 

Astron. Lett. 27 (2001) 799-808. 
[5] D. Markovic, "Bound near-equatorial orbits around 

neutron stars," arXiv:astro-ph/0009450. 
[6] T. Damour and G. Schafer, "Higher-order relativistic 

periastron advances and binary pulsars," Nuovo Cim. B 

101 (1988) 127-176. 
[7] G. Schafer and N. Wex, "Second post-Newtonian 

motion of compact binaries," Phys. Lett. A 174 (1993) 

196-205. 



[8] G. Schafer and N. Wex, "Erratum: Second 

post-Newtonian motion of compact binaries," Phys. 
Lett. A 177 (1993) 461(E). 
[9] G. Scharf, "Schwarzschild geodesies in terms of elliptic 
functions and the related red shift," arXiv : 1101 . 1207 
[astro-ph.GA] . 

[10] C. Konigsdorffer and A. Gopakumar, "Post-Newtonian 
accurate parametric solution to the dynamics of 
spinning compact binaries in eccentric orbits: The 
leading order spin-orbit interaction," Phys. Rev. D 71 
(2005) 024039, arXiv:gr-qc/0501011. 

[11] M. Tessmer, J. Hartung, and G. Schafer, "Aligned 

spins: Orbital elements, decaying orbits, and last stable 
circular orbit to high post-Newtonian orders," Class. 
Quant. Grav. 30 (2013) 015007, arXiv: 1207.6961 
[gr-qc] . 

[12] B. Carter, "Global structure of the kerr family of 

gravitational fields," Phys. Rev. 174 (1968) 1559-1571. 

[13] E. E. Flanagan and T. Hinderer, "Evolution of the 

Carter constant for inspirals into a black hole: Effect of 



5 



the black hole quadrupole," Phys. Rev. D 75 (2007) 

124007, arXiv : 0704 . 0389 [gr-qc] . 
[14] C. M. Will, "Carter-like constants of the motion in 

newtonian gravity and electrodynamics," Phys. Rev. 

Lett. 102 (2009) 061101. 
[15] G. W. Hill, "Motion of a system of material points 

under the action of gravitation," Astronom. J. 27 

(1913) no. 646-647, 171-182. 
[16] I. G. Izsak, "A note on perturbation theory," Astronom. 

J. 68 (1963) 559. 
[17] K. Aksnes, "On the Use of the Hill Variables in 

Artificial Satellite Theory," Astronomy and 

Astrophysics 17 (1972) 70. 
[18] R. H. Boyer and R. W. Lindquist, "Maximal Analytic 

Extension of the Kerr Metric," J. Math. Phys. 8 (1967) 

265-281. 

[19] A. Merloni et al, "On gravitomagnetic precession 
around black holes," Mon. Not. R. Astron. Soc. 304 
(1999) 155-159. 

[20] A. Gopakumar and G. Schafer, "Gravitational wave 
phasing for spinning compact binaries in inspiraling 
eccentric orbits," Phys. Rev. D 84 (2011) 124007. 



[21] M. Schneider and C. Cui, Theoreme iiber 
Bewegungsintegrale und ihre Anwendung in 
Bahntheorien, vol. 121 of A Theoretische Geoddsie. 
Verlag der Bayerischen Akademie der Wissenschaften, 
Miinchen, Germany, 2005. 

[22] A. Sommerfeld, Atombau und Spektrallinien, vol. 1. 
Friedr. Vieweg & Sohn, Braunschweig, 7 ed., 1951. 

[23] T. Damour, P. Jaranowski, and G. Schafer, "Dynamical 
invariants for general relativistic two-body systems at 
the third post-Newtonian approximation," Phys. Rev. D 
62 (2000) 044024, arXiv : gr-qc/9912092. 

[24] S. Gillessen et at, "Monitoring stellar orbits around the 
massive black hole in the galactic center," ApJ 692 
(2009) 1075. 

[25] Y. Kato et al, "Measuring spin of a supermassive black 
hole at the galactic centre implications for a unique 
spin," MNRAS Letters 403 (2010) L74-L78. 

[26] W. Schmidt, "Celestial mechanics in kerr spacetime," 
Class. Quant. Grav. 19 (2002) 2743. 

[27] R. Fujita and W. Hikida, "Analytical solutions of 
bound timelike geodesic orbits in kerr spacetime," 
Class. Quant. Grav. 26 (2009) 135002. 



